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ABSTRACT 

PSR B 1259-63 is a 48 ms pulsar in a highly eccentric 3.4 year orbit around the young massive star LS 2883. 
During the periastron passage the system displays transient non-thermal unpulsed emission from radio to very 
high energy gamma rays. It is one of the three galactic binary systems clearly detected at TeV energies, 
together with LS 5039 and LS I H-61 303. We observed PSR B 1259-63 after the 2007 periastron passage with 
the Australian Long Baseline Array at 2.3 GHz to trace the milliarcsecond (mas) structure of the source at three 
different epochs. We have discovered extended and variable radio structure. The peak of the radio emission is 
detected outside the binary system near periastron, at projected distances of 10-20 mas (25^5 AU assuming 
a distance of 2.3 kpc). The total extent of the emission is ^50 mas (~120 AU). This is the first observational 
evidence that non-accreting pulsars orbiting massive stars can produce variable extended radio emission at AU 
scales. Similar structures are also seen in LS 5039 and LS I -1-6 1 303, in which the nature of the compact object 
is unknown. The discovery presented here for the young non-accreting pulsar PSR B 1259-63 reinforces the 
link with these two sources and supports the presence of pulsars in these systems as well. A simple kinematical 
model considering only a spherical stellar wind can approximately trace the extended structures if the binary 
system orbit has a longitude of the ascending node of 57 ^ -40° and a magnetization parameter of cr 0.005. 
Subject headings: gamma rays: stars - pulsars: individual (PSR B 1259-63) - radio continuum: stars - stars: 
emission-Une, Be - stars: individual (LS 2883) - X-rays: binaries 
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1. INTRODUCTION 

The binary system PSR B1259-63/LS 2883' is formed by 
a young 48 ms radio pulsar in an eccen tric orbit of 3.4 years 
around a massive main-sequence star (iJohnston et al.l [l992l 
Il994i) . The parameters of the system are shown in Ta- 
ble [T] The spectral type of the massive star, 09.5 Ve, and 
some of the binary par ameters have been recently updated by 
iNegueruela et al.l (1201 li) . who obtained a distance to the sys- 
tem of 2.3 ± 0.4 kpc. Close to the periastron passage the sys- 
tem displays no n-thermal unpulsed e mission that has been de- 
tected in radio ( | Johnston et al.l2005'). X-rays (Comins ky et alJ 
11994 IUchivamaetal.1 12009: Chernvakova et al. 2009), hard 
X-rays up to 200 keV (iGrove et al.i il995i) . and very 
high energy (VHE: 0.1-100 TeV) 7-rays above 380 GeV 
dAharonian et 31.11200 5. 2009). In the range -0.1-100 GeV 
strict upper limits were ob tained by EG RET close to the 1994 
periastron passage (Tava ni et al.l 119961) . and the source has 
not yet been observed by AGILE and Fermi close to perias- 
tron. The VHE emission is variable on orbital timescales, 
and is interpreted as the result of inverse Compton upscatter- 
ing of stellar UV photons by relativistic electrons, which are 
accelerated in the shock between the relativistic wind of the 
young non-accreting pulsar and the w ind of the stella r com- 
panion (se e Maraschi & Treve s1ll9811: iTa vani & Aronsi (T997l 
iKirk et alJ[T999ciDubusii2006i:IBogovalov et alji2008i. and ref- 
erences therein). 

The high-mass binaries LS 5039 and LS I -1-6I 303 have also 

* Star 2883 in the catalog of Luminous Star's in the Southern Milky Way 
(Stephenson & Sanduleak|[l97lh . The use of SS 2883 should be avoided, see 
INegueruela et alj J201 III 



been detected at VHE, and show a broadband spectral energy 
distribution (SED) similai- to that of PSR B1259-63/LS 2883 
(iDubusli2006 ). In contrast to X-ray binaries, all three sources 
have the peak of the SED at MeV-GeV energies. For these 
reasons, they can be considered gamma-ray binaries. How- 
ever, the nature of the compact objects in LS 5039 and 
LS I -1-61 303 is unknown because their m asses are not 
well con strained by the system mass functions dCasares et alj 
2005a. bl: [Aragona et al.l 120091) . and no pulsations have been 
found. These systems have been extensively observed at 
VHE during several orbital cycles, while the observations of 
PSR B 1259-63, the only system with a confirmed pulsar, are 
scarce due to the long orbital period. Any observational link 
between the three gamma-ray binaries would shed light in the 
understanding of this kind of systems. 

The radio emission fro m the PSR B 1259-63 system is 
described most recently in IJohnston et al.l (|2005), which in- 
cludes multiwavelength data from ATCA observations ob- 
tained during the periastron passages of 1994, 1997, 2000, 
and 2004 (hereafter we will refer to the epoch of each pe- 
riastron passage as r). The emission has two non-thermal 
components. The first one is pulsed emission with a flux den- 
sity of —2-5 mJy at 2.5 GHz and a nearly flat spectral index, 
which disappears approximately from 16 days prior to perias- 
tron (t-16) to 15 days after periastron (t+15). This has been 
interpreted as an eclipse of the pulsar when it crosses behind 
the equatorial circumstellar disk present around the massive 
star (Melatos et al. 1995). The second component is transient 
unpulsed synchrotron emission that appears at t-20 days and 
shows two peaks centered around r-10 and r+20 days, with 
flux densities at 2.5 GHz up to — 15-20 mJy and —30-50 mJy, 
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TABLE 1 

Parameters of the Pulsar, the Massive Star and the Binary System. 



Parameter 


o 1 1 

Symbol 


Value 


Keterence 


Pulsar period 


P 


47.762506780(2) ms 


1 


Period derivative 


p 


2.276554(2) x 10"'^ 


1 


1^ n Qi'c*r*tf*n cf"i erf* 
V'lidJ. ax. LCI 1 A LIU uUC 




^ y 10^ vr 
/\ L\f y 1 


2 


Surface magnetic field 


B 


3.3X10 Cj 


L 


Spindown luminosity 


Esp 


8 X 10 " erg s ' 


3 


Spectral type 




09.5 Ve 


4 


Effective temperature 


Teff 


27500-34000 K 


4 


Surface gravity 


log g 


3.7^.1 


4 


Radius 


Ri 


8.1-9.7 Rq 


4 


Optical luminosity 


Lopt 


2.4 X 10 ° erg s ' 


4 


r^ass 




31 Mq 


4 


Distance 


d 


2.3 ± 0.4 kpc 


4 


Mass function 


mi) 


1.53 M0 


5 


Terminal wind velocity 


Voo 


1350 ±200 km s-i 


6 


Orbital period 


/'orb 


1236.72432(2) days 


1 


Reference epoch 


To 


MJD 48124.34911(9) 


1 


Semimajor axis 




7.2 ± 0.4 AU 


4a 


Orbit inclination 


i 


22°2±1.4 


4a 


Eccentricity 


e 


0.8698872(9) 


1 


Argument of periastron 




138.°6659(1) 


1 


Longitude of ascending node 


n 


-40° 


See the text 


Proper motion (rigth ascension) 


COS S 


—1.4 it 2.7 mas yr"' 


7 


Proper motion (declination) 




-3.2 ±1.9 mas yr"' 


7 



References. — (l)IWang"erani2004l; (2)IWex et al][T993 : (3)lManchester etlTI 
[1995; (4) Negueru ^et al.ll20m (5) IJohnston et aHfT^^l (6) IMcCoUuiiJfT^ ; (7) 
IZacharias et al. 200^ 

Note. — The values in parentheses refer to the uncertainty in the last digit at la- 
level. 

" Derived from (4). 



respectively. After the post-periastron peak, the flux density 
of the unpulsed emission decreases continuously, and it has 
been detected up to ^ r+lOO days. This transient emission 
remains optically thin during the outbursts. 

The broadband transient emission of PSR B 1259-63/ 
LS 2883 is produced around periastron, when strong inter- 
action is produced between the stellar and pulsar winds. The 
shocked material is contained by the stellar wind behind the 
pulsar, producing a nebula extending away from the stel- 
lar companion. Along this adiabatically expanding flow, the 
accelerated particles produce synchrotron emission from ra- 
dio to X-rays (Tavani & Arons 1997; Kirk et al. 1999; Dubus 
i2006t iTakata & Taamll2009D . The expected morphology de- 
pends on the magnetization parameter of the pulsar wind, a, 
defined as the upstream ratio of magnetic to kinetic energy. 
These models predict that the radio emission extends up to 
several AU (corresponding to several milliarcseconds, mas, 
at 2.3 kpc), and that its structure should be variable on or- 
bital timescales. This radio behavior has not been tested in 
PSR B 1259-63, which only emits over a period of a few 
months every 3.4 year. 

Here we present the first VLBI radio images of 
PSR B 1259-63, obtained close to its 2007 periastron passage. 
The high-resolution images at three different orbital phases 
provide a direct view of the small-scale morphology of the 
source, which is comparable to those previously observed in 
LS 5039andLSlH-61 303. 

2. OBSERVATIONS AND DATA REDUCTION 

PSR B 1259-63 was observed with the AustraUan Long 
Baseline Array (LB A) at 2.3 GHz (13 cm) on three epochs: 
2007 July 28 (run A), 2007 August 17 (run B), and 2008 



June 6 (run C). The LBA observations were performed with 
five antennas of the array: Parkes, ATCA, Mopra, Hobart (not 
present in run C) and Ceduna. The observational parameters 
of each of the ~ 10 hr runs are shown in Table |2] The small 
number of antennas provides a rather poor Mv-coverage that 
makes the data calibration and imaging more difficult than 
for arrays with more elements, because we have less base- 
line redundancy. The data were recorded at a bit rate of 
512 Mbps per telescope distributed in eight sub-bands (four 
for each right- and left-handed polarization) with a bandwidth 
of 16 MHz, each of them correlated using 64 frequency chan- 
nels, two-bit sampling, and 2 s of integration time. Hobart 
and Ceduna recorded at 256 Mbps (only two sub-bands for 
polarization). The data were correlated at Swinburne Univer- 
sity using the DiFX software corr elator (Deller etaLl 120071) 
without applying pulsar gating or binning. 

The observations were performed using phase referenc- 
ing on the calibrator J 1337-6509 (B 1334-649), which has 
an angular separation of 4.°0 from PSR B 1259-63 and 
was correlated at aj2ooo.o = 13*'37™52!'4443 and (5j2ooo.o = 
-65°09'24."900. This reference position from the global 
VLBI solution 2006d_astro^ has an uncertainty of 13 mas. 
The cycle time was 6 minutes, spending half of the time on the 
phase calibrator and the target source alternatively. The total 
flux density of the phase calibrator was 457 ± 11, 367 ± 20, 
and 491 ±6 mJy for runs A, B, and C, respectively. The 
source J0538-4405 (B0537-441) was used as a fringe finder 
for runs A and B, and J 1337-6509 (B 1349-439) was used for 
run C. No astrometric check source was observed during the 
runs. 

^ http://lacerta.gsfc.nasa.gov/vlbi/solutions/2006d/2006d_apr.src 
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Fig. 1. — LBA images of PSR B 1259— 63 at 2.3 GHz. North is up and east is to the left. The dates and the days after the periastron passage (r) are quoted at 
the top of each panel. The synthesized beam is displayed in the rectangle on the bottom-right comer of each image. The red crosses mark the region where the 
pulsar should be contained in each run (see the text). As a reference, the size of the major axis of the orbit of PSR B 1259— 63/LS 2883 is shown in the first panel. 
For each image, the displayed contours start at 3iT and increase by factors of 2'^^^; the la rms close to the source in each image from left to right is 0.30, 0.66, 
and 0.15 mJy beam"' . 



TABLE 2 

ObservationalParameters for Each Run. 



Run MJD On-source Time No. of Antennas r-Fdays" Orbital Phase" (°) Tnie Anomaly (°)'' 6(°)''' 



A 54309.04-54309.46 4.65 5 1.1(2)-1.5(2) 0.00087(14)-0.00120(14) 9.1(1.4)-12.5(1.4) 191 

B 54328.97-54329.40 3.75 5 21.0(2)-21.4(2) 0.01698(14)-0.01733(14) 98.3(3)-99.1(3) 279 

C 54623.29-54623.67 4.20 4 315.3(2)-315.7(2) 0.25496(14)-0.25527(14) 165.97(1)-165.99(1) 346 



Note. — ^ is the mean true anomaly plus 180° . 
" Periastron passage at r = MJD 54307.9710(1), from the fourth orbital solution obt 

The data reduction was performed in AIPS ^*. Total electron 
content models based on GPS data obtained from the CD- 
DIS data archive'* were used to correct phase variations due 
to the ionosphere. Standard instrumental corrections were ap- 
plied (parallactic angle, instrumental offsets, and slopes be- 
tween and within bands and bandpasses). Fringe fitting on the 
phase calibrator was performed with the AIPS task FRING, 
and the solutions were applied to the target source. A self- 
calibrated model of the calibrator was used as an input model 
for CALIB, which was used to reduce the amplitude instabil- 
ities on timescales greater than 10 minutes. The data were 
averaged in frequency and time, and clean images were pro- 
duced. For run B only, a single round of phase self-calibration 
was applied, to mitigate the residual atmospheric phase insta- 
bilities, which were more noticeable at this epoch. The fi- 
nal images were produced using a natural weighting scheme 
(robust 5 within the AIPS task IMAGR). For run B, robust 
2 and tapering were applied to avoid the presence of possi- 
ble unreliable high-resolution features due to sidelobes of the 
synthesized beam. Self-calibration slightly affects measured 
properties such as extent and position angle (PA.), while it 
preserves the morphology. 

3. RESULTS 

The resulting VLBI images at 2.3 GHz are shown in Fig- 
ure [U Extended emission is detected at distances up to 50- 
55 mas (120-130 ±20 AU at 2.3 ±0.4 kpc) during the two 
runs shortly after the periastron passage. The emission be- 

^ The NRAO Astronomical Image Processing System. 
|http://www. aips.nrao.edu/ 

^ The Cmstal Dynamics Data Information System http://cddis.nasa.gov/ 



in lWangetai] 120041) . 

comes gradually fainter from the peak toward the northwest, 
and no individual components have been found. The PA. of 
the extended emission with respect to the peak is ^ -67° 
for run A and ^ -50° for run B. The emission in run C, 
315 days after the periastron passage, is dominated by a point- 
like source of a few mJy. 

Like all VLBI arrays, the absolute flux calibration of the 
LBA relies on noise calibration injection and so the absolute 
flux values reported in Table |3]should be taken as uncertain at 
the ^ 10% level. The flux densities of runs A and B are com- 
patible with previous ATCA observations at the correspond- 
ing orbital phases (the ATCA data from the current observa- 
tions were not correlated as an independent array). The flux 
density in run C is compatible with the flux density of the pul- 
sar. This is expected considering the lack of unpulsed emis- 
sion at t+ 150 and r + 180 in previous ATCA observations 
jJohnston et al.l l2005'). 

The phase-referenced observations allow us to obtain rel- 
ative astrometry between runs. Since no astrometric check 
source was observed, we do not have real measurements of the 
astrometric uncertainties, and we will use the formal errors of 
a Gaussian fit obtained with JMFIT within AIPS. Given the 
relatively large calibrator-target separation of ^ 4°, we ex- 
pect an additional systematic component to the position error 
due to the unmodeled ionosphere of 1-5 mas (iBrisken et alj 
l2000h . As a reference position for the plots we use the peak 
position of run C, aj2ooo.o = 13''02"47?6435(1) and (Sj2ooo.o = 
-63 °50'08 "636(1), which we consider to represent the pul- 
sar position at MJD 54623.48. Assuming a mass of the neu- 
tron star of 1 .4 Mq and a stellar mass of 3 1 ±5 Mq (see Ta- 
ble [T] for the system parameters), the mass function provides 
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Fig. 2. — Computed trajectoiy of the nebular flow in the past. The contour plots are the same as in Figure [T](left and center). The green elhpse represents the 
counterclockwise orbit. The longitude of the ascending node, Q, is set to —40° . The different magnetizations used are displayed in the bottom-left panels. The 
axes units are AU (100 AU ~ 43 mas). 



TABLE 3 

Source Parameters for the Images Shown in Figure[T1 



Run 


(mJy) 


(mJy beam"' ) 


Shpbw 
(mas) 


P.A.HPBW 

(°) 




Size 


P.A. 

(°) 


Aa 
(mas) 


A<5 
(mas) 


Sepai'ation 


(mas) 


(AU) 


(mas) (AU) 


A 


19.9 ± 1.4 


10.4 ± 0.2 


28.9 X 26.1 


-49 


50 


120 ± 20 


-67 


11.3 ±0.4 


14.0 ± 0.5 


(10-20) ± 3 (24^6) ± 7 


B 


46.7 ± 1.0 


32.7 ± 0.4 


31.0 X 25.1 


-78 


55 


132 ±22 


-50 


4.2 ±0.1 


11.3±0.1 


(5-14) ±3 (12-31) ±7 


C 


3.0 ± 0.4 


2.8 ± 0.4 


50.3 X 25.1 


-14 


<2.8 


<6.7 ± 1.1 




0.0 ± 0.6 


0.0 ± 1.1 





Note. — The columns are mn label, total and peak flux density at 2.3 GHz, synthesized beam (HPBW) parameters, size and P.A. of the extended emission, the position offset 
of the peak of the emission from the reference position (position of mn C), and the range of possible separations between the peak and the pulsar position. 



/ = 22 .°2 and the semimajor axis of the pulsar orbit is 7.2 AU, 
or 3.1 mas for a distance to the system of 2.3 kpc. The red 
crosses in Figure [T] mark the region where the pulsar should 
be located in each epoch. Their centers are placed at the posi- 
tion of run C corrected for proper motion at the corresponding 
epochs (MJD 54309.25 for run A and MJD 54329.18 for run 
B). Since we do not know the orientation of the orbit in the 
sky, we plot as error bars the projected orbital separation of 
the pulsar with respect to run C. The error bars also include 
the 1(7 uncertainties on the mass of the star (and hence the 
distance uncertainty), on the astrometry of run C, and on the 
offset due to the proper motion. Finally, we have included the 
astrometry in runs A and B to compute the range of possible 
separations between the peak of the emission and the pulsar, 
as shown in Table [3] We consider all possible values of the 
longitude of the ascending node (fi), which determines the 
orientation of the orbit in the plane of the sky. 

4. KINEMATICAL INTERPRETATION 

The radio morphology at a given epoch depends on the spa- 
tial distribution of synchrotron emitting particles and their 
emission processes. Given the limitations of our data (only 
two images, and without accurate astrometry), we have used a 
simple kinematical model to check if it can trace the extended 
structures detected. 

We have considered the shock between the relativistic pul- 
sar wind and a spherical stellar wind. The shock is produced 
at the standoff distance, the region where the pulsar and stellar 



wind pressures balance, as described in iDubusI (120061) . The 
evolution of the nebular fl ow after the shock is described in 
iKennel & Coronitil (119841) . This should be considered as a 
first approximation to the much more complex hydrodynamic 
behavi or of a shocked flow in a binary system, as shown in 
Bogov aiov et alj ( 120081) ■ In the Kennel & Coroniti approxi- 
mation of a non-turbulent adiabatically expanding flow, the 
flow speed depends only on the magnetization parameter cr 
when assuming ct 1. This allows us to compute the past 
trajectory of the flow produced behind the standoff distance, 
which depends on the components separation along the orbit, 
the mass lo ss rate of the st ar (0.6 x 10"^ Mq yr"' for a typi- 
cal 09 star: IVink et al.ll200 0). the terminal wind velocity, Voo, 
and the spin-down luminosity of the pulsar, E^p. With these 
restrictions, the only free parameters are the longitude of the 
ascending node, il, which describes the orientation of the or- 
bit, and the magnetization parameter, cr. Projecting the past 
trajectory of the flow on our VLBI images of runs A and B, 
we found that the best match with the detected morphologies 
is obtained for ~ -40° and a magnetization parameter of 
cr ~ 0.005, assuming an orbital inclination of / = 22.°2. The 
obtained trajectories are shown in Figure |2]for three different 
values of cr. We note that the uncertainty in the distance to 
the system scales the size of the contours, but not the orbit 
and the trajectories, which are computed in AU. The range of 
magnetizations in the plots has been chosen to approximately 
show the effect of keeping cr = 0.005 and changing either the 
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distance from 1.9 to 2.7 kpc, either from -35 to -45°, or a 
variation of the product M Voo of two orders of magnitude. 

The parameters of the circumstellar equatorial disk of the 
star are uncertain, but considering M^, = 5 x 10"** M yr~', 
dJohnston et al.|[T996h . and a wind velocity of ^10 km s"', 
the product Me ^oo.e is ^100 times smaller than the spher- 
ical wind contribution. If this equatorial component domi- 
nates during certain orbital phases, the flow trajectory would 
be closer to the orbit than the dashed models in Figure |2]for 
some specific regions of the more recent part of the trajec- 
tory. However, we do not attempt to accurately model this 
circumstellar disk, as its density, velocity, and crossing time 
are unknown. 

We note that the astrometric errors can be of the order of 
^10 AU (see above). We also emphasize that this simple 
model cannot account for the complex magnetohydrodynam- 
ical turbulences of the real flow and should be considered as a 
first approximation to constrain a. Previous studies have as- 
sumed a magnetization parameter a of around 0.01-0.02 for 
this pulsar (see Tava ni & Aronsll99 7: Dubus 2006), consider- 
ably higher than our best-fit value, although comparable lower 
values for a have been suggested fo r other systems (e.g., the 
Crab pulsar: iKennel & Coroniti|[T984i) . 

5. DISCUSSION AND CONCLUSIONS 

The results presented in this Letter show that the particle ac- 
celerator within PSR B1259-63/LS 2883 can produce a flow 
of particles emitting synchrotron radiation that can travel sev- 
eral AU. The total projected extent of the nebula is ^ 50 mas, 
or 120 ± 20 AU, and the peak of the emission is clearly dis- 
placed from the binary system orbit (see Figure [T] and Ta- 
ble [3). Similar morphologies and displacements have been 
found in the other two known gamma-ray binaries, LS 5039 
and LS I -h 61 303, although for smaller sizes and on shorter 
timescales (Ribo et al.' 2008; Moldon et al. '2 010l;lM"assi et alj 
|2b04; Dhawan et al. 2006; Albert et al. 2008). 

There is an ongoing debate on the nature of the com- 
pact object and particle acceleration mechanisms in LS 5039 
and LS I H-61 303. Although initial ly suggested to be 
accreting/ejecting microquasar systems (iParedes et al.ll200^ 
iMassi et ail 120041) . they are now thought to contai n young 
non-accreting pulsars d Maraschi & Treves 1981; Dubus 
iDhawan et all 12006 ; Ribo et ak 2008; Moldon et al 



which can explain their multiwavelength emis 



sion ("Sierpowska-Bartosik & Torre sl2008 : iCe r utti et al.l2008l ; 
[Bogovalov et al. 2008; Zdziarski et al.l 20101) . However, 
there are three basic issues that are not well understood 
for LS 5039 and LS I H-61 303, and comparisons with 
PSR B1259-63/LS 2883 can help to clarify the situation. 
First, the putative pulsar properties of these two sources are 
unknown, as no pulsations have been detected for these sys- 
tems. The lack of pulsations can be explained by the in- 
tense stellar wind that produces an extremely high absorp- 
tion for these small orbits, a 3.9 day orbit with separations 
between 0.1 and 0.2 AU for LS 5039, and 2 6.5 days with 
separations of 0.1-0.7 AU for LS 1 H-61 303 (iCasares et al.l 



l2005allbl ; lAragona et al.ll2009h . As a reference, the separation 
for PSR B1259-63 is in the range 0.9-13.4 AU, and the pul- 
sations disappear for distances below ^1.6 AU. Second, it 
is not clea r if the massive ste llar wind can confine the pul- 
sar wind (iRomero et al.ll20()7b . VLBl images, like the ones 
presented here, can shed light on the shock conditions and 
geometry. Third, the observed SED and variabihty at GeV 
energies is not well understood (lAbdo et al.ll2009allbt iTorresI 
120101) . Fermi and AGILE are observing PSR B 1259-63 for 
the first time during the 2010 periastron passage, and are pro- 
viding GeV data that can be compared with those already ob- 
tained for LS 5039 and LS 1 +61 303. In this context, the 
high-resolution VLBl radio observations presented here es- 
tablish a common link to test the similarities between the three 
systems. 

In conclusion, our results provide the first observational ev- 
idence that non-accreting pulsars orbiting massive stars can 
produce variable extended radio emission at AU scales. Sim- 
ilar structures are also seen in LS 5039 and LS 1 +61 303, 
in which the nature of the compact object is unknown be- 
cause the detection of pulsations is challenging. The dis- 
covery presented here for the young non-accreting pulsar 
PSR B 1259-63 reinforces the link with these two sources 
and supports the presence of pulsars in these systems as well. 
Planned LBA observations of PSR B 1 259-63 during the 20 1 
periastron passage will allow us to provide a full comparison 
with the behavior observed in these sources, which have been 
extensively monitored during several orbital cycles. We have 
also shown that the orientation of the orbit and the magnetiza- 
tion of the pulsar can be inferred from VLBl observations of 
the source. Several images at different orbital phases cover- 
ing a wider range of true anomalies will allow for a complete 
modeling of the orbital changes of the extended emission. 
Finally, accurate VLBl observations of the pulsed emission 
during several orbits can provide the pulsar trajectory, from 
which we can directly obtain the proper motion of the binary 
system, the inclination and fl of the orbit, and the distance to 
the system. 
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